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INTRODUCTION 

In  the  past  dynamic  behavior  of  infinite  periodically  laminated  medium  has  been  stu¬ 
died  extensively.  A  review  of  the  literature  on  exact  and  approximate  analyses  of  this  prob¬ 
lem  can  be  found  in  [1,2]. 

—  In  this  paper  we  use  the  stiffness  method  that  was  presented  in  [2]  to  study  disper¬ 
sion  of  waves  in  a  laminated  plate.  In  this  approach  each  lamina  is  divided  into  several 
sublayers  and  the  displacement  distribution  through  the  thickness  of  each  sublayer  is 
approximated  by  polynomial  interpolation  functions  in  such  a  way  that  displacements  and 
tractions  are  continuous  across  the  interfaces  between  adjacent  sublayers.  Details  of  the 
method  can  be  found  in  {2,3].  Here  we  summarize  the  pertinent  equations  and  discuss 
numerical  results  obtained  for  particular  systems. 

This  study  is  motivated  by  a  desire  to  model  wave  propagation  in  a  continuous  fiber- 
reinforced  laminated  plate.  If  the  wavelength  is  long  compared  to  the  fiber  diameters  and 
spacing,  then  each  lamina  can  be  modeled  as  a  homogeneous  transversely  isotropic 
medium  with  the  symmetry  axis  parallel  to  the  fibers.  The  overall  effective  elastic  proper¬ 
ties  of  such  a  medium  can  be  calculated  from  the  fiber  and  matrix  properties  by  using  an 
effective  modulus  theory  [4,5].  Such  an  assumption  has  been  made  in  this  paper.  Thus 
the  dispersion  characteristics  of  guided  waves  in  a  layered  anisotropic  plate  has  been  stu¬ 
died  here.  '  .  n  fj  • , 
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GOVERNING  EQUATIONS 

We  consider  time  harmonic  waves  in  a  plate  composed  of  m  laminae.  For  simplicity 
in  analysis  it  will  be  assumed  that  each  lamina  is  transversely  isotropic  with  the  symmetry 
axis  aligned  with  either  the  x-  or  the  y-axis  (Fig.  1).  This  assumption  is  not  necessary  for 


r.v  u  — 


the  development  of  the  equations,  but  it  is  made  here  to  keep  the  algebra  as  simple  as  pos¬ 
sible  for  the  anisotropic  problem  at  hand.  Under  this  assumption  the  wave  propagation 
problem  reduces  to  two  uncoupled  ones:  plane  strain  in  which  the  displacement  com¬ 
ponents  are  ux,  0,  u„  and  SH  or  antiplane  strain  when  the  only  nonzero  displacement 
component  is  uy.  So  in  the  following  we  shall  treat  these  two  separately. 


Plane  Strain 

Consider  the  ith  lamina  bounded  by  z=Z;_i  and  z=zi+1.  The  stress  strain  relation  in 
this  lamina  will  be  given  by 
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where  a y  and  are  the  stress  and  strain  components,  respectively,  and  we  have  written 

'X„=2ca.  Note  that  if  y-axis  is  the  axis  of  symmetry,  then  c$=c$  and  (c fp-clj*). 


The  problem  is  then  equivalent  to  that  of  an  isotropic  one. 

In  order  to  get  good  numerical  results  each  lamina  will  be  divided  into  several  sub¬ 
layers,  Ms,  say.  Within  the  jth  sublayer  we  will  choose  a  local  coordinate  with  the  origin  at 
the  mid-plane  and  x}>  yjt  zj(  parallel  to  the  global  x,y,z  axes,  respectively.  Let  2^  be  the 
thickness  of  this  sublayer.  Denoting  u^  to  be  the  displacement  at  a  point  in  the  jth  lam¬ 


ina  we  write 
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where  4,  (n=l,...,4)  are  cubic  polynomials  in  the  local  coordinate  Zj  given  by 
h  =  ~(2-3r;1+fjj3),  f,  = 


f3  =  f<  =  -^(-l-'?j+'?/+'7jS) 


(3) 


Here  r)^—zjh  ■  and  Uj,  uj^,  Xj,  ai  are  the  values  of  ux,  u„  an,  and  £r„  at  the  jth  node.  These 
nodal  values  of  the  displacement  and  traction  components  are  functions  of  Xj  (-  x)  and  t. 

In  this  paper  it  will  be  assumed  that  the  time  dependence  is  of  the  form  e-iwt,  u>  being  the 
circular  frequency.  The  factor  e-*"*  will  be  dropped  in  the  sequel. 

The  equations  governing  the  nodal  generalized  coordinates  {uj,  Xj,  Wj,  <7j}  are  obtained 
using  Hamilton’s  principle.  For  details  see  [3].  It  is  found  that  the  amplitudes  of  the  gen¬ 
eralized  nodal  variables  satisfy  the  matrix  equation 


(k<(K  j  -  iks[K  1  -  ks(K  ]  +  ik|K  ]  +  |K  ])  {Q  >  =  0 
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where  we  have  written 


(4] 


<Q}  *  {Q  }eikx  (5) 

{Q}  being  the  vector  of  all  the  generalized  coordinates.  Matrices  [K  ],  etc.,  have  been 

defined  in  [3],  Note  that  [K  ],...,(k  ]  depend  linearly  on  ui1. 
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For  nontrivial  solution  the  determinant  of  the  coefficient  matrix  of  {Q0}  must  be 
zero.  This  is  the  dispersion  equation  to  solve  for  k  for  given  uj. 

Antiplane  Strain 

The  derivation  given  above  is  for  the  plane  strain.  The  antiplane  strain  case  can  be 
considered  in  a  similar  manner  and  that  results  in  a  much  simpler  frequency /wave  number 
equation.  In  this  case  we  assume 


Here  r=c$du  ^)/5z.  The  corresponding  eigenvalue  problem  for  the  wave  number  is 
(-k*[Kj  +  [K  ])  {Q  }  =  0 

As  illustrative  examples  we  solved  equations  (4)  and  (7)  for  a  fiber-reinforced  (0 
homogeneous  plate  and  a  three-layered  (0  *  /90  *  /0  * )  plate  with  or  without  interface 
layers.  These  results  are  discussed  in  the  next  section. 


(6) 

(7) 

) 


NUMERICAL  RESULTS  AND  DISCUSSION 

As  an  application  of  the  technique  described  above,  we  considered  a  fiber-reinforced 
plate  when  the  fibers  are  aligned  along  the  x-axis  (0  * ).  The  properties  of  the  plate  are 
given  in  Table  1.  As  mentioned  before,  for  propagation  in  the  0*(90*)  direction  the  prob¬ 
lem  is  tractable  analytically.  Figure  2  shows  the  frequency  spectrum  for  propagation  in  the 
x-direction.  Also  shown  in  this  figure  is  the  conipariron  with  the  exact  solution.  It  is 


Fig.  5.  First  three  real  branches  of  the  frequency  spectrum  for  a  0  *  /90  *  /0  *  plate  with 
soft  interface  thin  layers. 


Finally,  Fig.  6  shows  the  same  effect  on  the  first  three  branches  of  the  frequency  spec¬ 
trum  of  SH  waves. 
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Fig.  6.  First  three  real  branches  of  the  frequency  spectrum  of  SH  waves  for  a  0  *  /90  *  /0  * 
plate  with  soft  interface  thin  layers. 

CONCLUSION 

Dispersion  curves  are  presented  for  propagation  of  in-plane  and  out-of-plane  (SH) 
waves  propagating  along  the  fibers  in  a  fiber-reinforced  plate.  These  curves  differ  signifi¬ 
cantly  from  those  in  an  isotropic  plate.  It  is  also  shown  that  the  dispersion  curves  for 
waves  propagating  in  a  cross-ply  (0  *  /90  *  /0  * )  laminate  with  interface  layers  differ  signifi¬ 
cantly  from  those  for  a  laminate  without  interface  layers  and  both  of  these  are  different 
than  the  ones  for  a  single-layered  fiber-reinforced  plate.  The  most  significant  features  of 
the  curves  for  a  plate  with  interface  layers  are  the  lowering  of  the  cut-off  frequencies  and 
their  slopes.  This  should  have  significant  effect  on  the  overall  plate  response  and  would 
allow  ultrasonic  characterization  of  interface  properties. 
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seen  that  the  comparison  is  excellent.  In  this  figure  n  = - -  — —  -  -  and  7  =  - .  It 

2^n/(c ti/p)  o-  2* 

is  seen  that  as  k— *00  the  slopes  of  the  first  symmetric  and  antisymmetric  modes  tend  to 
the  ratio  >/C55/p,  where  V  R  is  the  Rayleigh  wave  velocity  in  the  x-direction.  For  the 
stiffness  method  we  used  15  sublayers. 

Table  1.  Properties  of  0  *  and  90  *  laminae,  and  the  interface  layer.  All  the  stiffness  are 
in  the  units  of  10nN/m  *. 
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Fig.  2.  Dispersion  curves  for  the  real  and  imaginary  wave  numbers.  Solid  lines  are 

obtained  by  the  present  method.  Symbols  x  (symmetric)  and  o  (antisymmetric) 
are  the  exact  calculations. 


The  next  problem  considered  as  a  sandwich  plate  with  0  *  /90  *  /0  *  configuration. 
Figure  3  shows  the  dispersion  curves.  These  differ  appreciably  from  the  ones  shown  in  Fig. 
2.  It  is  seen  that  the  cut-off  frequencies  are  lowered  as  well  as  the  slopes  of  the  curves. 
Branches  of  nonpropagating  modes  are  also  quite  different. 


Figure  4  shows  the  dispersion  curves  for  a  0  *  /90  *  /0  *  plate  with  two  interface  layers 
between  the  two  outer  laminae  (0* )  and  one  inner  lamina  (90* ).  It  is  assumed  that 
thickness  ratio  of  each  interface  layer  and  each  adjacent  lamina  is  0.1.  Cut-off  frequencies 
are  seen  to  be  lower  than  those  in  Fig.  3.  To  see  the  effect  of  interface  layers  on  the  first 
three  branches  of  the  frequency  spectrum  we  plotted  renormalized  H  vs.  7  in  Fig.  5.  In 

this  figure  O  =  - - —  •  -  111  and  7  =  — .  It  is  interesting  to  note  the  slowing  down 


of  the  waves  as  the  wavelength  becomes  of  the  order  of  the  thickness  of  the  interface 
layers. 


Fig.  3.  Dispersion  curves  for  0  *  /90  *  /0  *  plate.  Here  H  is  the  total  thickness  of  the 

plate,  fl  =  - — -—===,  and  7  = 

2*V(cw/?)  cr 
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